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DETERMINATION  OF  VISCOEIASTIC  MODEL  CONSTANTS  FROM 
DYNAMIC  MECHANICAL  mOPERTIES  OF  LINEAR  VISCOELASTIC  MATERIALS 

ABSTRACT 

A  semi-analytical  method  of  determining  the  generalized  Voigt  model  which 
represents  the  dynamic  mechanical  properties  of  a  linear  viscoelastic  material 
over  a  range  of  frequencies  of  three  decades  is  presented.  This  model 
representation  is  shown  to  be  equivalent  to  the  differential  operator  formulation 
of  the  linear  viscoelastic  stress-strain  law.  The  method  is  applied  to  complex 
creep  compliance  data  for  N.B,  S. poly isobutylene  at  22  different  temperatures. 

In  general,  the  compliances  calculated  from  the  models  differ  from  the 
experimental  data  by  less  than  5^- 

A  summeiry  of  spring- das hpot  model  theory  is  presented  in  Appendix  A.  The 
equivalence  of  the  differential  operator  stress-strain  relation  and  the 
generalized  Voigt  model  is  demonstrated  in  Appendix  B. 


5 


Page  intentionally  blank 


Page  intentionally  blank 


Page  intentionally  blank 


TABLE  OF  COWTEWTS 


ABSTEACT  . 

TABLE  OF  CONTEMTS . 

INTRODUCTION  . 

LINEAR  ISOTHERMAL  VISCOELASTIC  STRESS- STRAIN  RELATIONS  .  . 

MODEL  REPRESENTATION  . 

COMPLEX  CREEP  COMPLIANCE . 

DETERiyCrNATION  OF  MODEL  CONSTANTS  FROM  DYNAMIC  CREEP 

COMPLIANCE  DATA . 

APPLICATION  OF  PROCEDURE  TO  N.B.S.  POLYISOBUTYLENE  .  .  .  . 

SUMMARY  AND  CONCLUSIONS . 

ACKNOWLEDGMENTS . 

REFERENCES  . 

APPENDIX  A 

SUMMARY  OF  SPRING-DASHPOT  THEORY . 

APPENDIX  B 

EQUIVALENCE  OP  DIFFERENTIAL  OPERATOR  0  -  e  RELATION  AND 
THE  GENERALIZED  VOIGT  MODEL . 

DISTRIBUTION  LIST . 


Page 

3 

5 

7 

7 

8 

9 

11 

13 

'Ti'Z. 

JJ 

33 

34 

37 


41 

43 


5 


Page  intentionally  blank 


Page  intentionally  blank 


Page  intentionally  blank 


IMTRODUCTION 


A  niomber  of  commonly  used  materials  such  as  plastics,  rubbers,  fibers, 

etc.  exhibit  viscoelastic  properties  such  as  creep  and  relaxation.  A  case 

in  point  is  the  solid  propellant  rocket  fuel  where  elastic  analyses  are 

inadequate  for  relatively  long  time  conditions,  such  as  slump  during  storage, 

for  extremely  short-time  dynamic  conditions,  such  as  impact  and  sudden 

acceleration.  Problems  such  as  these  are  currently  being  solved  using  techniques 

(l)  (2)* 

of  the  theory  of  linear  viscoelasticity.'  However,  the  application  of  these 

techniques  is  dependent  upon  the  formvilation  of  the  linear  viscoelastic  stress- 
strain  laws  from  experimental  data  in  such  a  manner  as  to  represent  the  mechanical 
properties  of  the  actual  material. 

Presented  in  this  report  is  a  semi-analytical  method  of  determining  the 

generalized  Voigt  model  which  will  accurately  represent  the  dynamic  mechanical 

properties  -  of  a  linear  viscoelastic  substance  over  a  range  of  frequencies  of  three 
(3) 

decades.  This  model  representation  is  equivalent  to  the  differential  operator 
formulation  of  the  linear  viscoelastic  stress-strain  laws. 

The  method  employed  is  based  upon  the  analogy  between  linear  electrical 

(M  (5) 

networks  and  spring-dashpot  mechanical  systems.'  '  a  transfer  function 

is  constructed  to  represent  the  generalized  Voigt  model.  The  asymptotic 

approximation  method  originally  devised  by  Baum  is  used  in  conjunction  with 

a  least  squares  computer  program  to  determine  the  viscoelastic  model 

constants  from  dynamic  creep  compliance  data.  This  method  can  be  applied 

to  viscoelastic  materials  exhibiting  a  wide  range  of  dynamic  mechanical 

behavior  and  represents  a  significant  improvement  over  previous  model  fitting 
(3) 

techniques.  ' 

LINEAR  ISOTHERMA.L  VISCOELASTIC  STRESS- STRAIN  RELATIONS 

In  the  theory  of  linear  isothermal  viscoelasticity,  a  viscoelastic  solid 
is  specified  by  the  existence  of  a  relation  connecting  stress,  a,  and  strain,  e: 

p[c]-«[ej  (1) 


* 

Superscripts  indicate  references  at  the  end  of  this  paper. 
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where  P  and  Q  are  linear  differential  operators  with  respect  to  time.  Thus, 


"  =  Po  +  Pi  dt  "  Pp  712 

U.0 


p. 


^  dt^ 


(2) 


Q  = 

yj 


+ 


d 

dt 


c 


+ 


+ 


(5) 


The  p's  and  q's  are  constants  which  represent  the  mecheinical  properties 
of  the  material  at  a  given  temperature.  These  constants  cannot  be  meastired 
directly,  but  must  be  calculated  from  experimental  data  in  which  both  the 
stress  and  strain  are  known  as  functions  of  time.  Although  the  stress-strain, 
o  -  €,  relations  for  a  linearly  viscoelastic  material  may  also  be  represented 
by  hereditary  integrals,  it  has  been  foxmd  that  the  determination 

of  the  stresses  in  a  viscoelastic  body  subjected  to  transient  loading  is 
considerably  simplified  if  Eqs.  (l),  (2)  and  (5)  are  used  to  specify  the 
material  behavior. 


MODEL  REPRESENTATION 

The  expression  P  = 
representation  of  a  mechanical  model,  consisting  of  a  network  of  springs  and 
dashpots.  A  summary  of  the  theory  of  spring-dashpot  models  is  presented  in 
Appendix  A.  The  model  used  mechanical  properties  of  materials 

in  this  report  is  the  generalized  Voigt  model. 


[' 


e 


is  often  visualized  as  a  mathematical 
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FIG.  i.  THE  GENERALIZED  VOIGT 

MODEL 
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This  model  exhibits  both  instantaneous  elasticity  and  long-term  creep  and 
is  mathematically  better  suited  to  describe  experimental  dynamic  creep  data 
than  is  the  generalized  Maxwell  model.  The  operational  form  of  the  stress- 
strain,  a  -  e  relation,  for  this  model  is  then  given  by 


e(t)  = 


q(t) 


n 


n..  — 

■u  ut; 


_1  \  * _ 1 

^  ^  E,  +  1 


i=l  i 


'i  dt- 


a(t) 


rriMPT.TSY  OPWIPP  POMPT.TAMrin 


The  steady  state  response  of  this  model  to  a  sinusoidal  stress  is  given 


Viv 


(5) 


lent 


where  the  time  dependence  of  both  e  and  a  is  of  the  form  e  This  ratio,  J*, 

is  known  as  the  complex  compliance  with  J*  and  j"  called  the  storage  and  loss 
compliances  respectively.  The  non-zero  Imaginary  function,  J"(a)),  implies  a 
phase  difference  between  the  oscillations  of  stress  and  strain.  Both  J*  and  J" 
functions  of  cd  which  can  be  measured  experimentally  over  the  frequency  range  of 
interest.  Substitution  of  steady  state  sinusoidal  variation  in  Eq.  (4), 
which  means  replacement  of  the  operator  with  io),  and  comparing  with  Eq.  (5) 


yields 


n 


J*(a))  =  ^  + 


0 


Tlo^ 


list] 


1-1 

O”-*- 


and  therefore 


E 


Ti 

”  ■  E^ 


n 


„  2  2  2 
E-  +  n. 

'1  “ 


„  2  2  2 
Hi  +  n  /n 


(6) 


/  rrN 


n 


•n 


n 


E. 

J 


E, 


'0 


/  .  „  -d  d  d 

E^  +  oi  n  . 
J-J-  J  'J 


are 
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1  "^1^ 
J"(£d)  =  -i-  +  . 


2^ 
<« 


+  *  •  •  + 


2  2  2 
E  +  n  a> 
n  'n 


(8) 


li 


"o'" 


M 


The  complex  compliance  may  also  be  written  as  the  ratio  of  two  polynomials. 
In  Eqs.  (2)  and  (5),  let  r  =  m;  this  corresponds  to  a  viscoelastic  material 
exhibiting  instantaneous  elasticity.  Then  choose  Pq  =  1  and  =  0;  this  allows 
for  the  presence  of  an  isolated  dashpot  which  permits  long  term  viscous  flow. 

Then  if  p  =  e^'" 

e  =  J*(aj)e^‘^^ 

and  from  (2)  and  (5)  we  have 


(9) 


ml  iajfc 


and 


but 


F  LCfJ=  +  P2(i“)  +  •••  +  Pjj^(M“Je 

Q  r€l=  +  qg(la)^  + 


icot 


(10) 

(11) 


Therefore 


0^(10)) 


1  +  p^ico  +  P2(ia))“  +  ...  +  Pjjj^(ia3)" 


'^r 


1(1)  4.  n 


r%  I  -f  ...  1 

qg\ia>; 


f  •  •  • 


/ .  \  m 
q^(lcu; 


(12) 


Now  let  IcD  =  s.  Then 


H(8)  =  ^  = 


1  +  p^s  + 


2  m 

ns  +  . . .  +  p  s 


;=:  m 

q^s  +  q^s  +  . . .  +  qjjjS 


(15) 


(7) 

Cauer, '  '  Alfrey ,  Gross  and  others  have  shown  the  analogy  between  electrical 
networks  and  mechanical  models.  From  this  analogy,  H(s)  is  defined  as  the 
system  transfer  function,  and  exhibits  the  same  properties  as  the  transfer 
function  for  an  RC  network.  The  equivalence  of  the  differential  operator 
o  -  €  representation  and  the  generalized  Voigt  model  is  presented  in 
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Appendix  B  along  with  a  discussion  of  the  properties  of  the  system  transfer 
function  h(s).  Another  convenient  form  of  the  transfer  function  H(s)  can 
be  derived  by  letting  ...  X,^_i  be  the  roots  of  the  equation 


.  _2  .  ,  .  .m-1 

To+t-nncr  II  II  .  -  .  n  Khp  r'nnf.R  of*  't'.hR  pnna+.inn 

- - - - -  - - 

2  111 

1  +  p^s  +  PgS  +  . . .  +  p^s  =  0 


/  -1  I.  \ 


(15) 


The  properties  or  H(s),  as  founa  in  Appenaix  a,  require  that  the  and 
be  negative  real  constants.  Then  H(s)  may  be  written  in  the  form 


T./  \ 

ni^s; 


TJL  /  ^  \ 

j-l.s;  = 


k^s  -  ii^Hs  - 


u  -  d  ; 

m 


!(S  -  X^)(s  -  x^)  •••  (s  - 


where  k  is  a  positive  real  constant.  This  form  of  H^s)  is  useful  in 
calculating  p(s)  and  q(s)  once  the  X^.  and  |i^.  have  been  determined. 


/  -1  /'\ 
Vio; 


DETERMINATION  OF  MDDEL  CONSTANTS  PROM 
DYNAMIC  CREEP  COMPLIANCE  DATA 


Having  shown  that  the  generalized  Voigt  model  is  equivalent  to  the 
differential  operator  description  of  viscoelastic  material  behavior,  the 

c»T\-r\v*r\-r»v*T  o+ca  'trolnoa  n-P  n  flnrl  W  P»r»P  r'#am T?  "hn  •Pn+'.  mian+.n  1  /^a+*o  rpVica 

V  »<VX  WCV-.  kj  V./X  ij  .  4.J  ,  w  .k  ww  j.  -i.  w  w.n.^Wdk  >,/«..».  X  •  XXXW 

J  J 

dynamic  creep  compliance  data  consists  of  values  of  J’  and  tP'  at  different 
values  of  cd.  From  the  relations  i 


|j*l  =  +  (J")^  (17) 

and 

^  -1  J" 

0  =  tan  -jY" 

the  values  of  1  the  absolute  value  of  the  complex  compliance,  and  0,  the 
phase  angle,  are  calculated. 

/n\ 

The  graphical,  semi-analytical  method  reported  by  Welch, and 
originally  applied  to  circuit  synthesis  is  then  used  to  determine  the 
approximate  location  of  the  poles  of  the  transfer  function,  H(s).  This 
method  consists  of  (l)  construction  of  a  graph  of  0  versus  log  co  and  (2) 
construction  of  a  graph  of  log  | J*|  vs  log  co  from  the  experimental  data. 
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Careflil  study  of  the  phase  angle  plot,  along  with  some  knowledge  of  the  effect 

of  various  smcings  of  poles  and  zeros  on  the  phase  angle  as  a  function  of 

(9) 

frequency,  results  in  an  initial  placement  of  the  poles  and  zeros.  The 
possible  location  of  these  poles  and  zeros  is  then  verified  on  the  grajdi  |  J*j 
vs  0).  These  poles  are  located  at  s  =  Xj  in  the  complex  s  plane.  These  values  of 


'j 


are  then  substituted  into  Eqs.  (7)  and  (8)  along  with  the  relation 


E./n  (Appendix  B)  to  yield 
J  J 


(18) 


J*(a>) 


/n  1 


0 


+  U) 


^2^  ^2 


+  O) 


X  /n 
n  'n 

K  n 

X  +0) 
n 


+ 

Now  let  =  A;  =  Bj 
Then  Eq.  (15)  becomes 
J'  (oj)  =  A  - 


'/nr 


'/n 


.  2  2 

Xj^  +0) 


4" 


n 


+  CD 


,  2  2 
X  +  CD 
n 


-1 


=  C 


(19) 


(20) 


— J 
-r  CD 


X  C 
n  n 

^n 


(21) 


J^'(cd)  =  -  + 

'  ''  (D 


cdC, 


cdC 


n 


n 


CD 


Eqs.  (21)  are  then  solved  for  the  values  of  A,  B,  and  Cj  using  experimental 

data  for  J'  and  J"  at  from  8  to  I6  different  values  of  cd.  The  solution  of 

this  overdetermined  system  is  obtained  using  a  least  squares  program  on  the 

BEL  ORDVAC  digital  computer.  This  program  applies  Legendre's  principle,  as 

outlined  in  Reference  10,  to  solve  the  normal  equations  derived  from  the 

overdetermined  system  of  linear  equations.  The  computer  uses  the  values  of 

X.  and  calculates  positive  A,  B,  and  C  along  with  the  values  of  J'  and  J"  at 
J  J 

the  different  values  of  cd  for  comparison  with  the  input  data.  Since  the  relative 

experimental  error  is  generally  constant,  the  least  squares  program  minimizes 

(ill* 

the  relative  error  in  the  fitting  procedure . '  ' 


This  procedure  has  proven  more  effective  than  the  method  employed  in  Reference  11 
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study  of  these  results  provides  indications  as  to  possible  motion  of 

the  poles  and  zeros  of  H(s)  to  improve  the  data  fit.  For  example,  if  the 

first  choice  of  the  poles  X.  is  in  error,  the  constants  C.,  as  calculated 

J  J 

by  least  squares,  are  negative.  These  negative  Cj  do  not  correspond  to  a 

spring  dashpot  model.  Possible  error  in  only  one  =  X^^  results  in  only 

C,  being  negative,  thus  providing  a  clue  to  the  source  of  the  trouble.  Once 
K 

the  A,  B,  and  Cj  are  known,  the  values  of  and  can  be  determined  from 

Eqs.  (l8)  and  (20).  The  problem  of  generating  p(3)  and  q(s)  from  the  Ej  and 

T) .  is  most  easily  approached  through  the  equivalence  of  two  forms  of  H(s). 

J 

Let  im  =  s  in  Eq.  (6)  and  substitute  Eq.  (20)  to  get 


J*(s)  =  H(s) 


Equating  Eqs.  (22)  and  (l6),  we  find 


and  the  are  the  roots  of  the  equation. 


(22) 


(25) 


Solution  of  Eq.  (25)  using  the  BRL  ORDVAC  digital  computer  by  means  of  successive 

subdivision  of  the  Intervals  between  X .  and  X . . n ,  provides  the  values  of  the  . 

j  J+x  ^i 

APPLICATION  OF  PROCEDURE  TO  N.B.S.  POLYISOBUTYLENE 

(12) 

Published  data  for  N.B.S.  polyisobutylene,'  '  Table  No.  I,  are  used  to 
provide  an  example  of  the  application  of  this  procedure.  The  data  for  22 
different  temperatures  are  fitted,  since  this  represents  a  wide  range  of 
viscoelastic  behavior.  The  results  of  these  calculations  are  tabulated  in 
Tables  No.  II  and  No.  Ill  and  are  Illustrated  in  Figures  10  thru  l8.  Figures  2 
thru  9  show  the  graphs  of  0  vs  log  cs  and  |  J* |  vs  log  a>  for  four  characteristic 
temperatures;  -40.4°C,  -29.8°C,  -9.9°C,  and  +50.0°C.  These  0  vs  log  oo  graphs 
indicate  considerably  different  viscoelastic  behavior,  and  therefore  the 
location  of  the  Xj,  as  well  as  their  spacing,  requires  some  knowledge  of  pole 
and  zero  interaction. 
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TABLE  I  SUMMARY  OF  DYNAMIC  MECHANICAL  DATA  ON  POLY  ISO  BUT  YLENE  * 


T*C 

Freqoeocy*  cps 

30 

40 

45 

60 

72 

80 

100 

IW 

200 

210 

2W 

400 

600 

800 

1000 

1400 

1500 

2000 

2800 

3000 

4000 

4200 

5100 

-44.6 

<!iii>/<lyne(X10^) 

0.151 

0.145 

0.136 

0.123 

0.132 

0.119 

0.111 

0.108 

0.0M5 

0.345 

0.244 

OJOO 

0.166 

0.142 

0.118 

A  no.u 

Q.OQl 

0.069? 

0.0551 

0.0506 

0.0368 

n  mno 

-40.4 

cm«/dyiie(X10-») 

0.308 

04M 

6.2i2 

6.186 

oJ«‘ 

0.167 

0.149 

0.137 

0.129 

oim' 

0.823 

0.617 

0.400 

0.411 

0436 

0484 

0430 

0.177 

0.151 

0.129 

0.107 

0.0021 

0.0651 

0.0601 

-34.: 

j.  cni»/dynt{X10^) 

0.848 

0.702 

0.571 

0.478 

0.387 

0.308 

0.254 

0.248 

0.216 

0.188 

0.165 

0.139 

2.20 

1.59 

1.17 

0.056 

0405 

0.660 

0.566 

0.467 

0.367 

0.304 

0.268 

0414 

0.182 

0.130 

0.107 

-29.8 

cmVdyneCXIO-*) 

1.41 

1.19 

0.097 

0.852 

0.656 

0.512 

0.440 

0.398 

0468 

0478 

0420 

0.153 

3.96 

2.92 

2.07 

1.70 

1.43 

1.13 

1.02 

0.796 

0.639 

0.553 

0.486 

0.415 

0.336 

0448 

0.226 

0J19 

-25.0 

jw  cm’ /dyne  (X 10^) 

2.54 

2.17 

1.81 

1.48 

1.25 

l.Oi 

0.K25 

0.703 

0.625 

0.526 

0.422 

0.320 

0.302 

*.* 

6.79 

4.98 

3.66 

2.97 

240 

1.80 

1.58 

1.26 

1.02 

0.880 

O.TH 

0.622 

0.558 

0.421 

0.350 

-19.9 

j»cmVdyne(X»0-*) 

6  56 

5.23 

3.S3 

3.42 

2.83 

2.25 

1.88 

1.52 

1.24 

1.03 

0.947 

0.792 

0.646 

0428 

OJSg 

... 

11.8 

8.91 

6.56 

548 

4.20 

3.14 

2.63 

2.08 

1.63 

1.37 

1.23 

0.977 

0.893 

0.687 

0.634 

-14.8 

jLemVdynetXIO-*) 

11.5 

8.96 

7.48 

€.17 

5.55 

4.51 

3.64 

2.99 

2.38 

1.85 

1.47 

... 

0.949 

0.752 

20.1 

14.8 

12.2 

104 

8.76 

6.03 

5.47 

4.42 

3.49 

2.65 

1.91 

1.40 

1.10 

-  9.9 

yi,cni»/dyiie(XIO-») 

19.7 

14.8 

12.0 

10.0 

8.70 

.  7.06 

5.87 

4.71 

3.55 

2.78 

2.43 

2.18 

1.62 

1.25 
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On  the  graphs  of  log  J*  vs  log  o)  are  shown  the  straight-line 
approximations  made  hy  the  method  of  asymptotic  slopes^ using  the 
calculated  poles  and  zeros.  These  straight-line  approximations  do  not  completely 
follow  the  calculated  curve  in  some  cases;  this  effect  results  from  the 
interaction  of  errors.  The  approximation  error  at  an  isolated  pair  of  critical 

/o\ 

values  is  given  by^ 

2. 

;)  (24) 


A\) 


rw) 


X 

-  log  /2  -  i  log  (1  +  -n  , 


4. 

It  a  pair  of  critical  values  in  which  the  ratio  —  is  fairly  large  is  followed 

by  a  pair  in  which  it  is  relatively  smaller,  the  two  errors  may  interact.  However, 

such  fluctimtions  are  of  little  importance,  since  in  most  cases  a  knowledge  of 

the  phase  angle  will  permit  the  approximation  to  be  made  with  sufficient  accuracy 

for  the  use  of  the  computer  program  mentioned  previovisly.  In  90^^  of  the  trials 

made  in  fitting  the  polyisobutylene  data,  the  first  choice  of  the  X.  provided  a 

J 

fair  fit  for  the  experimental  data.  These  models  generally  were  within  +5^  in 
the  middle  range  of  data  with  some  having  errors  of  +10^  at  the  extremes.  Slight 
adjustment  of  the  resulted  in  models  which  fit  the  data  at  92^  of  the  data 
points  with  an  accuracy  of  +5^  and  with  only  one  data  point  deviating  by  more  than 
10%. 

Of  significant  importance  is  the  fact  that  once  the  X.  have  been  determined 

for  a  particular  set  of  experimental  data,  introduction  of  new  data  which 

are  more  or  less  consistent  with  the  old  does  not  require  selection  of  new  X., 

J 


but  only  requires  the  use  of  the  computer  program  to  recalcvilate  the  other  model 
parameters.  As  an  example,  recent  data  on  polyisobutylene  at  T  =  which 

differ  from  the  Fitzgerald  data  by  approximately  10^  were  fitted  using  the  same 
Xj.  Then  by  taking  the  geometric  mean  of  the  two  sets  of  model  constants,  the 
model  which  best  fits  all  the  experimental  data  was  determined  and  the  resultant 
values  of  J*  and  are  displayed  in  Figures  (19)  and  (20). 
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SUmiABY  AND  CONCLUSIONS 


The  results  show  that  the  combination  of  the  semi- analytical  and  least 


V/  A 


•f’-japhni  nni=».c« 


r6a,dily  produces  models  'bhu'b  fit 


experimental  data  over  a  3  decade  range  of  frequency.  The  models  generated 
by  this  method  can  describe  a  wide  range  of  viscoelastic  behavior  with  a 
nominal  accuracy  of  +5^.  The  application  of  this  method,  however,  is  limited 
to  materials  that  behave  in  such  a  manner  as  to  have  |  J*(a))  j  a  monotonically 
decreasing  function  of  cjo  whose  slope  is  always  greater  than  minus  one  when 


plotted  on  a  log- log  graph.  If  the  material  does  not  quite  meet  this 
criterion,  as  in  polyisobutylene  at  T  =  -44.6°C,  +80.0°C,  and  +99‘9°C,  then 
the  input  data  to  the  computer  program  must  be  adjusted  to  meet  this  requirement. 
Otherwise,  the  least  squares  program  will  yield  negative  Cj  in  attempting  to 
fit  the  data.  After  making  the  necessary  adjustments,  the  resultant  model  will 
still  closely  approximate  the  material  behavior. 


As  more  and  more  experimental  data  become  available,  the  models  generated 
by  this  method  can  easily  be  extended  to  fit  a  wider  range  of  frequencies. 
Since  these  models  describe  the  physical  behavior  of  viscoelastic  materials 
with  an  accuracy  that  is  of  the  order  of  experimental  errors,  they  can  be 
used  to  formulate  the  stress-strain  laws  necessary  to  the  solution  of 
analytical  problems.  Laboratory  experiments  can  then  be  conducted  using  the 
actual  materials,  and  a  comparison  of  analytical  resiilts  with  experimental 
results  can  then  be  made. 
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APPENDIX  A 


SUMMARY  OP  SEHING-nASHPOT  MODEL  THEJORY 

Combinations  of  spring  and  dashpot  elements  which  form  the  spring- 
dashpot  models,  have  been  found  to  provide  a  convenient  means  for  constructing 
a  differential  operator  in  order  to  approxi^te  observed  viscoelastic  behavior. 
The  basic  elements  of  these  models,  as  well  as  the  two  basic  combinations  of 
these  elements  are  shown  in  Figure  1. 


ELEMENTS 


The  spring  (Fig.  l.a)  represents  the  elastic  response  of  the  Baterial,  and 
if  a  and  e  represent  the  stress  and  strain,  then 

a(t)  =  Ee(t)  (l) 

where  E  is  a  constant  called  the  modulus  of  the  spring. 

The  dashpot,  consisting  of  a  piston,  cylinder,  and  Newtonian  fluid, 
represents  the  viscous  flow  of  the  material.  The  stress-strain  relation  is 

<j(t)  =  n  ^  (2) 

where  t)  is  a  constant  representing  the  viscosity  of  the  fluid  in  the  dashpot. 

These  two  elements  may  be  connected  in  two  ways:  (l)  in  series  (Fig.  l.c) 
Isiown  as  the  Maxwell  model,  and  (2)  in  parallel  (Fig.  l.d)  known  as  the  Voigt 
or  Kelvin  model. 

The  stress-strain  relation  for  the  Maxwell  model  is 

de  _  1  ^  0 
It  E  dt  Tj 
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and  for  the  Voigt  model  is 

a(t)  =  Ee(t)  +  n  ||  (^) 

Comparison  of  the  response  of  the  Maxwell  model  and  the  Voigt  model  with 
actual  materials  has  shown  that  nKsdels  consisting  of  more  elements  eire  required 
to  accurately  descrihe  the  mechanical  hehavior  of  viscoelastic  materials  over  a 
significant  range  of  frequencies. 

GENERALIZED  MAXWELL  MODEL 

This  model  is  an  extension  of  the  simple  Maxwell  model  which  can  he  tised  to 
describe  a  viscoelastic  material. 


I 


FIG.  2.  THE  GENERALIZED  MAXWELL 

MODEL 


The  stress-strain  relations  derived  from  Figure  2  are  given  hy 


for 


6=e„  =e  =e„  +€ 

®0  ’’u  'll 

i  =  1,  2,  . . .n 

0  =  a„  +  a  +0T+0„  +  ...+ct 

Eq  no  1  2  n 


=  nr 


de 

dt 


^  ^  '"JL  _d  _l~j 

dt  [^E^  dt  m  *^1 


(5) 
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_  ri  +  -ii  n 

dt  Eg  dt  ^2J  ^ 


^  =  ri  _§  +  -lift 

*  dt  E  dt  Ti  n 

n  'n 

Combination  of  these  equations  yields  the  operational  form; 


GENERALIZED  VOIGT  MODEL 


(6) 


This  model  is  an  extension  of  the  simple  Voigt  model,  and  is  the  model 
used  to  describe  the  mechanical  properties  of  the  materials  considered  in 
this  report.  This  model  exhibits  both  instanteineous  elasticity  and  long  term 


r-AArn  r-W^ 

■  *}i~*~*  I' . . 

^0  %  HH  HH 

Vi  Vz  'In 

FIG.  3.  THE  GENERALIZED  VOIGT 

MODEL 


creep  and  is  mathematically  better  suited  to  describe  experimental  dynamic 
creep  data  than  is  the  generalized  Maxwell  model. 

Tiie  stress-strain  relations  derived  from  Fig.  3  are  given  by; 

0 

de 

«  -  « 


a  =  ^ 


1  dt 


+  E,€- 
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a  =  E„€  +  x\ 

n  n  ' 


de 
_ n 

n  dt 


€=6  +e  +€,+€„+•••+€ 
Tip  1  2  n 


0  ''0 

Combination  of  Eqs,  (7)  yields. 


■(t)  =  . 


n 


z 


°  L  ^0  dl  ii  ®i  ^i  dt. 


a(t) 


(8) 
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AffHSWDiX  B 


B^mVALEMCE  OF  BIFFEREKTIAL  OPERATOR  a  -  €  RELATION 
AND  THE  GENERALIZED  VOIGT  MODEL 

The  use  of  the  generalized  Voigt  model  to  approximate  viscoelastic 
behavior  is  based  upon  the  equivalence  of  the  model’ s  resultant  differential 
operator  with  the  relation 

P[53=Q[ej.  (1) 

The  differential  operator  form,  for  the  generalized  Voigt  model,  as  found  in 
Appendix  A,  is  given  by  __ 

n 

6  (t)  -  ^  +  Y — j-  U(t).  (2) 

‘‘o  Ln„  4  +  n,  ^  J 

—  'o  uu  i=x  1  1  ai:  “• 

For  the  case  of  steady  state  sinusoidal  stresses  and  strains,  we  have  found 
that  for  the  model 

=  —  +  — —  4-  \  - i- -  .  (^\ 

Eg  ^  Ej  + 


Letting  j  =  ics 


j.(,)  =  + -i-  y.-.-i... 

Ej,  Z^E  +n,< 


We  have  also  shown  that  for  steady  state  sinusoidal  stresses  and  strains,  the 

-V'* 

differential  operator,  cj  -  e,  viscoelastic  relatldli  yields 

1  +  p,  its  +  p„(iu))^  +  ...  +  p  (io))”^ 

j^Cio.)  - - i - (5) 

q^ioD  +  q^Cics)  +  ...  +  q^(ia))‘“ 


T  ^4-  A  ^  ~  rrTU£.y» 

AWW  XCU  XXlCli 


J*(b)  =  H(s)  = 


P(6)  Z  ^1°  ^g°  •••  ^  V 


<lj_8  +  <lsS  +  .  .  .  + 


where  H(s)  Is  defined  as  thss  systeast  transfer  function. 


To  show  the  equivalence  of  Eqs.  (6)  and  (4),  let  X^,  ... 


be  roots  of  the  equation 


•  •  • 


The  roots  of  this  equation  must  be  distinct,  since  we  are  considering 
only  non-degenerate  models.  Then,  the  partial  fraction  expansion  for  H(8) 
is  given  by; 

m-1 


h(b)  =  _  +  ^  J 


q^s 


j=l 


(8) 


Now  let  m  -  1  =  n,  then 


n 


,/  1  ‘m  -t  \  1 

fu  -i-  ^  ^ 


Comparison  of  Eqs.  (9)  and  (4)  yields 


E  = 

0  P, 


m 


^0 


il> 


and 


p’  (l  w 

p(ijT 


which  then  yields 


E, 


(9) 


(10) 


(11) 


These  equations  relate  the  values  of  the  model  constants  to  the  material 
constants  in  the  stress-strain  law.  For  a  real  mechanical  system,  E,  and  q 
are  real  positive  constants  (or  zero).  Therefore,  the  are  real  and 
negative.  This  is  a  property  of  the  transfer  function,  H(s).  As  a  direct 
resialt  of  the  electrical  network  analogy,  the  following  are 

the  properties  of  the  system  transfer  function,  H(s)  for  the  generalized  Voigt 
model: 


(1)  all  zeros  and  poles  of  H(s)  occur  on  the  negative  axis  of  reals. 

(2)  all  zeros  and  poles  are  simple,  and  are  interlaced  in  the  pattern., 
pole-zero-pole-zero-. . . 

(3)  the  highest  critical  frequency  of  H(b)  is  a  zero;  H(oo)  is  a 
positive  real  number. 

(4)  the  lowest  critical  frequency  of  H(s)  is  a  pole  which  occurs  at 
the  origin,  s  =  0;  H(0)  =  oo . 
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